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Abstract — Based  on  energy  conservation  considerations 
we  study  the  nonlinear  dynamic  behavior  of  a quantum 
mesoscopic  circuit,  which  is  characterized  through 
parameters  of  inductance  and  capacitance.  Nonlinearity  is 
given  by  the  initial  conditions  of  magnetic  flux  and 


OC 

discreteness  charge  which  oscillate  in  the  interval  [-  0 , + 

],  being  the  magnetic  flux  normalized  by^ 

This  LC  circuit  with  quantized  electric  charge  is  excited 
by  energy  battery  that  can  produce  an  electrical 
discreteness  charge  on  the  capacitor.  The  dynamics  of  the 
mesoscopic  circuit  is  highly  nonlinear.  Our  results  show  for 
the  magnetic  flux  a nearly  square  wave  with  an  elongated 
period  when  compared  with  the  linear  case  and  a train  of 
narrow  pulses  for  the  discrete  charge. 
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I.  INTRODUCTION 

In  a series  of  articles  Li  and  Chen  [1,2]  and  Flores  et  al  [3, 
4,  5,  6,  7],  have  developed  a theory  of  quantum  electrical 
systems,  based  on  a treating  such  systems  as  quantum  LC 
circuits;  that  is,  electrical  systems  described  by  two 
parameters:  an  inductance  L,  and  a capacitance  C.  Such 
quantum  theory  of  circuits  is  expected  to  apply  when  the 
transport  dimension  becomes  comparable  with  the  charge 
carrier  coherence  length,  taking  into  account  both  the 
quantum  mechanical  properties  of  the  electron  system,  and 
also  the  discrete  nature  of  electric  charge.  Now,  we  propose 
a semiclassical  theory  of  quantum  electrical  circuits,  to 
obtain  useful  predictions  of  the  theory  from  very  simple 
calculations  of  energy  consideration  , and  to  push  the  circuit 
analogy  one  step  further,  generalizing  the  Heisenberg 
equations  of  motion  . 

The  semiclassical  theory  of  quantum  LC  circuits  [1]  starts 
from  the  quantum  Hamiltonian  of  the  LC  circuit  [1-7,  8,  9]. 
The  resulting  equations  become 
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which  we  can  put  as 

no: 

a+  C0q  sin(a)  = 0 (3') 

where  CO2  = 1 / LC  , OC  = Cf)  / cf)0 . The  equations  above  are 

considered,  mathematically,  as  classical  equations,  but  they 
include  quantum  effects,  the  quantized  nature  of  electric 

charge  through  the  parameter  fl  / qe  = (|)0  . These  equations 
are  highly  nonlinear,  but  they  reduce  to  the  usual  equations 
of  the  LC  circuit,  in  the  discrete  charge  limit,  qe  — » 0 . in 

this  paper  we  are  interested  in  the  highly  nonlinear  behavior 
of  the  LC  circuit  excited  by  a battery  source  with  initial 

conditions  (|)0  ^ 0 and  d((|)  / (|)0  ) / dt  ^ 0 , and 


We  will  resolve  the  nonlinear  equation  through  energy 
considerations  because  our  starting  point  is  the  Hamiltonian 
equation  (1) 


II.  DIFFERENTIAL  EQUATION  FROM 
ENERGY  CONSIDERATIONS 

At  the  initial  instance,  the  energy  of  the  mesoscopic  circuit 
is  the  sum  of  the  kinetic  and  potential  energy,  from  (4)  we 
obtain  : 

7 12 

E0  = CCo  + Ep  sin2(a0 12)  (4) 
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with  E0  = Eb=0  + Eb  , Ep  = 4C0q  =4/  LC  which  has 

O' 

been  normalized  and  OC  = Cf)  / (f)0  OC  = d((|)  / (f>0  ) / dt , 


fl  / qe  = (f)0  . Because  of  the  assumed  lossless  system  , the 

initial  energy  is  preserved  for  all  instances: 

'I. 

E0  = a2  (t)  + Ep  sin2  (a(t)  / 2) 

(5) 

Here,  Ep  / 2 = 2(jl)q  is  the  maximum  possible  (normalized) 

potential  energy  of  the  quantum  circuit,  being  attained  when 

it 

a = 0. 

Equation  (5)  refers  to  the  conservation  of  the  total 
mechanical  energy,  to  integrate  and  apply  the  initial 

n 

conditions  given  OC0 , OC0  , we  obtain 

da  o 

( — )2  = 2(Oofcos(a(t))-cos(a0(t))l  + ao  (6) 
dt 

The  analogy  between  a quantum  LC  circuit  with  discrete 
load  and  nonlinear  pendulum  is  straightforward  because 
both  systems  are  modeled  by  the  same  equations  (3'),  (6). 
These  equations  are  similar  to  that  obtained  for  nonlinear 
pendulum,  The  simple  gravity  pendulum  is  a famous  case 
study  in  classical  mechanics  that  leads  to  a nonlinear 
differential  equation  of  second  order  like,  equation  (3').  A 
solution  of  the  differential  equation  (6)  is  based  on  Jacobi 
elliptic  integrals  has  been  well  known.  There  exists  a great 
number  of  papers  and  textbooks  dealing  with  pendulums, 
the  reader  is  referred  to  [10], 

So  following  [10-13]  and  after  several  steps  algebraic  we 
can  obtain  the  normalized  magnetic  flux,  in  equation  (6) 
(see  also  appendix  A of  [14]),  where  we  have  incorporated 

the  initial  condition  of  da  / dt  . 


a(t)  = A = 2sin  1 
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, k = sin  (-±),  a0  = (<|>/  <f>0)|t=0 


(8) 


In  particular,  if  the  angular  velocity  is  zero,  (Xo=  0,  then 
Uq  = 1 therefore  equation  (22)  reproduces  the  known 
formula  for  the  nonlinear  pendulum,  [12,  13]. 

a(t)  = 2 sin-1  [Vksn(F(sin_1  (1),  k)  ± co0t,  k 

(9) 

71 

However  F(—  ,k)  = K(k)  is  the  first  class  elliptical 
integral  so 

a(t)  = 2 sin-1  [ Vksn(K(k)  ± co0t,  k) 

(10) 

Generally  the  first  equation  of  (8)  cannot  have  any  value,  it 
must  be  limited  because  for  a LC  circuit  the  amplitude 

ae  [-a0,a0]  so  0<k<l,and 


Uok  = k + ^T  = keff 
4con 


(11) 


which  also  keff  must  be  keff  < 1 . If  we  put  keff  ~ 1 , then 


go 

4w2 
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The  period  T is  the  time  required  to  complete  a cycle,  in  this 
case,  the  oscillation  period  T is  four  times  the  time  taken 

from  a = 0, (v  = 0) , a = a0, (v  = 1)  so 
co0/co  = 4K(k) 

(13) 

1 

the  parameter  k is  related  to  ao  and  CD0  as 

1 1 2 

k = a0 

2co0 

(14) 


III.  NUMERICAL  RESULTS 

Here  we  show  the  result  of  calculation  of  CO  as  a function 

of  a0  . For  simplicity,  we  assume  that  CO0  = 1 . For  small 
values  of  a = Cf)  / (|)0 , the  nonlinear  equation  (10)  becomes 

equal  to  a resonant  circuit  where  the  waveform  is 
completely  sinusoidal. 
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In  figure  1 , we  show  OC  = (|)  / (|)0  versus  t / T for 

OC0  = 71  / 2 . having  as  parameter 

1 f§. 

OC0  = (([))  / (f)0  (t  = 0)  = —0.5,  dashed  point-line,  0,  solid 
line  and  0,5  dashed  line  respectively.  The  magnetic  flux  is 
normalized  by  (j)0  = fl  / qe  , that  is  OC  = (f)  / (J>0 . 


U . fl 

having  as  parameter  OC0  = ((())  / (f)0  (t  = 0)  = —0.5,  dashed 
point-line,  0,  solid  line  and  0,5  dashed  line  respectively. 


At  small  oscillations,  the  energy  E0,  equation  (4)  of  the 
circuit  is  small  compared  to  the  maximum  possible 
potential  energy  Ep  = 4-CDq  . This  leads  to  a modulus  close 

to  zero,  for  which  the  Jacobi  elliptic  functions  can  be 
replaced  by  trigonometric  functions 

The  period  of  oscillation  of  the  mesoscopic  circuit  is 
constant  and  independent  of  the  initial  angular 

displacement  for  values  of  OC^—7ll  2 , and 

(dOC  I dt)  0=0  as  shown  in  Figure  1,  solid  line.  As  the 

initial  magnetic  flux  increases  (and  its  initial  time 
derivative  ),  the  oscillation  period  increases  considerably 

when  OC0  ~ 7T . That  is,  the  oscillation  becomes 

disharmonious.  For  this  reason,  once  the  movement 
becomes  dissonant  or  anharmonic  significantly,  the  period 
(frequency)  no  longer  remains  constant  but  lengthens 

fl 

increasingly  (shortens)  with  increasing  amplitude  of  OC0  as 

is  shown  in  figures  2 and  3.  Figure  3 is  an  enlarged  version 

□ 

of  Figure  2 near  the  value  a0  = 2 corresponding  to 
equation  (23)  with  k = 1 . 
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Fig.2:  Plot  of  CD  vs  OC0  for  (D0  = 1.  The  angular  frequency 

n 

reduces  to  zero  at  OC0  —2. 


CO 


n U 

Fig. 3:  Plot  of  CD  vs  0CQ  for  CD0  = 1,  1.99  < OC0  < 2 . The 

n 

angular  frequency  reduces  to  zero  at  OC0  = 2 . 


Fig.4:  The  amplitude  OC  and  the  angular  velocity  0C  are 
also  periodic  functions  for  0CQ  =0.99571,  : EB  =0.5.  At 

microwave  frequency  t / T = CD0t  X 271 X 1 0 8 , 

CX(t)  = (()/ ([)0  is  nearly  a square  waveform. 
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Figure  4,  shows  the  amplitude  OC  and  the  angular  velocity 

□ 

OC . They  are  also  periodic  functions  for  OC0  = 0.99571 , 


with  EB  = 0.5  , but  OC  is  almost  a square  wave  andOC  is  a 
train  of  narrow  pulses.  Here,  the  dynamics  of  the 

mesoscopic  circuit  is  highly  nonlinear.  Here,  OC  is 
proportional  to  the  discreteness  charge  Q of  the  quantum 
circuit 


Q 

c 


1 


= -0C 


(24) 

Waveforms  of  dOC  / dt  are  appropriate  for  design  of  bits 
sequence  for  digital  systems  at  high  frequency. 


IV.  CONCLUSIONS 

Based  on  energy  conservation  considerations  we  study  the 
nonlinear  dynamic  behavior  of  a quantum  mesoscopic 
circuit,  which  is  characterized  through  parameters  of 
inductance  and  capacitance.  Nonlinearity  is  given  by  the 
initial  conditions  of  magnetic  flux  and  discreteness  charge 
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which  oscillate  in  the  interval  [-OC0,  +OC0],  being  ao  the 


magnetic  flux  normalized  by  ft  / qe  = (|)0  and  it  is  close  to 

±71.  This  LC  circuit  with  quantized  electric  charge  is 
excited  by  energy  battery  that  can  produce  an  electrical 
discreteness  charge  in  LC  in  the  form  of  narrow  pulses.  The 
results  presented  also  intend  to  bring  the  student  of  physics 
and  engineering,  the  introduction  of  elliptic  integrals  and 
motivate  the  search  for  new  alternatives  in  the  design  of 
electronic  circuits  besides  to  solve  physical  and  applied 
mathematical  problems. 
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